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It is necessary to use algorithms to analyze gene expression 
data from DNA microarrays, such as in clustering and ma- 
chine learning. Previously, we developed the knowledge- 
based fuzzy adaptive resonance theory (KB- Fuzzy ART), a 
clustering algorithm suitable for analyzing gene expression 
data, to find clues for identifying gene networks. Leaf primor- 
dia form around the shoot apical meristem (SAM), which 
consists of indeterminate stem cells. Upon initiation of leaf 
development, adaxial-abaxial patterning is crucial for lateral 
expansion, via cellular proliferation, and the formation of flat 
symmetric leaves. Many regulatory genes that specify such 
patterning have been identified. Analysis by the 
KB-FuzzyART and subsequent molecular and genetic ana- 
lyses previously showed that ASYMMETRIC LEAVES1 (AST) 
and AS2 repress the expression of some abaxial-determinant 
genes, such as AUXIN RESPONSE FACTOR3 (ARF3)/ETTIN 
(ETT) and ARF4, which are responsible for defects in leaf 
adaxial-abaxial polarity in asl and asl. In the present 
study, genetic analysis revealed that ARF3/ETT and ARF4 
were regulated by modifier genes, BOBBER1 (BOB 7) and 
ELONCATA3 (£L03), together with AS1-AS2. We analyzed 
expression arrays with asl elo3 and asl bob1, and extracted 
genes downstream of ARF 3 /ETT by using KB-FuzzyART and 
molecular analyses. The results showed that expression of 
Kip-related protein (KRP) (for inhibitors of cyclin-dependent 
protein kinases) and Isopentenyltransferase (IPT) (for 



biosynthesis of cytokinin) genes were controlled by AS1- 
AS2 through ARF3/ETT and ARF4 functions, which suggests 
that the AS1-AS2-ETT pathway plays a critical role in con- 
trolling the cell division cycle and the biosynthesis of cyto- 
kinin around SAM to stabilize leaf development in 
Arabidopsis thaliana. 

Keywords: ASYMMETRIC LEAVES2 (AS2) • AUXIN RESPONSE 
FACTOR3/ETTIN • CDK inhibitors • Cytokinin • Shoot apical 
meristem. 

Abbreviations: CDK, cyclin-dependent kinase; CKI, cyclin- 
dependent kinase inhibitor; GCOS, Gene Chip Operating 
Software; KB-FuzzyART, knowledge-based fuzzy adaptive 
resonance theory; qRT-PCR, quantitative real-time reverse 
transcription-PCR; SAM, shoot apical meristem; ta-siRNA, 
trans-acting short interfering RNA; TSA, trichostatin A. 



Introduction 




Over the past several years, rapid advances in DNA microarray 
technologies, including a tiling array procedure, have made it 
possible to measure the expression levels of thousands of genes 
simultaneously and under different conditions (Yoshimura 
et al. 2011). These volumes of gene expression data can be 
used to reveal the mechanisms that drive life. To analyze this 
large quantity of data, however, requires the use of algorithms, 
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as used in clustering and machine learning. We previously de- 
veloped knowledge-based fuzzy adaptive resonance theory 
(KB-FuzzyART), which facilitates clustering analysis on the 
basis of a home-made list including a limited number of 
genes of interest, and applied the analysis of gene expression 
data to identify genetic networks that might be responsible for 
a physiological process (Takahashi et al. 2008). 

In general, the objective of a statistical or informatical ana- 
lysis is the enrichment of important information. In particular, 
we have selected the KB-FuzzyART clustering method in the 
present study. This method can enable us not only to enrich 
important information, but also to separate information into 
that which is known or unknown. We updated our previously 
used gene list (Gene-list-1) of Arabidopsis thaliana to make a 
new gene list (Gene-list-2). Using Gene-list-2, we applied 
KB-FuzzyART to two series of microarray data, to identify the 
genes downstream of ASYMMETRIC LEAVES1 (AS7) and 
ASYMMETRIC LEAVES! (AS2) of A. thaliana, both of which 
are involved in leaf development (Byrne et al. 2000, Iwakawa 
et al. 2002, Semiarti et al. 2001). 

Leaf primordia are formed at the periphery of the shoot 
apical meristem (SAM), which contains aggregates of indeter- 
minate stem cells. Homeostasis of SAM was recently shown to 
be controlled by multiple receptor protein kinases as well as the 
SHOOT-M ERISTEM LESS gene (Takano et al. 2010, Betsuyaku 
et al. 2011, Uchida et al. 2013). As each leaf grows, its morph- 
ology becomes established along three axes, the proximal- 
distal, adaxial-abaxial and medial-lateral axes. Adaxial-abaxial 
patterning at the initial stage, occurring in regions adjacent to 
the SAM, is critical for the lateral expansion of the lamina along 
the medial-lateral axis for formation of flat symmetric leaves 
(Waites and Hudson 1995, Byrne et al. 2000, Semiarti et al. 2001, 
Tsukaya 2006, Iwakawa et al. 2007, Bowman and Floyd 2008, 
Szakonyi et al. 2010, Moon and Hake 2011, Yamaguchi et al. 
2012). 

In A. thaliana, many genes for putative transcription factors 
have been identified as regulators of leaf development. Three 
members of the family of class III homeodomain-leucine zipper 
genes, such as PHABULOSA (PHB), specify the adaxial identity of 
lateral organs (McConnell and Barton 1998, McConnell et al. 
2001, Emery et al. 2003, Bao et al. 2004, Mallory et al. 2004). 
Members of the KANADI (KAN) and FILAMENTOUS FLOWER 
(FIL)/YABBY (YAB) gene families have been identified as abaxial 
determinants (Bowman and Smyth 1999, Eshed et al. 1999, 
Sawa et al. 1999, Siegfried et al. 1999, Eshed et al. 2001, 
Kerstetter et al. 2001, Bowman and Floyd 2008, Goldshmidt 
et al. 2008, Sarojam et al. 2010). The other components deter- 
mining abaxial cell identity include the AUXIN RESPONSE 
FACTOR3 (ARF3)/ETTIN (ETT) and AUXIN RESPONSE 
FACTOR4 (ARF4) genes (Sessions and Zambryski 1995, Kim 
et al. 1997, Sessions et al. 1997, Ulmasov et al. 1997, Pekker 
et al. 2005, Wu et al. 2008, Kelley et al. 2012). The results of 
investigations of these previously identified transcription fac- 
tors support the hypothesis that the specification of adaxial- 
abaxial polarity might be tightly coupled with lateral expansion, 



because inhibition of the development of either the adaxial or 
abaxial domain causes the appearance of filamentous leaves 
with the surface characteristics of either domain (Waites and 
Hudson 1995). Recently, it has been reported that succinic 
semialdehyde dehydrogenase is involved in leaf patterning 
(Toyokura et al. 2011). AUXIN RESISTANT3 also has a role in 
the formation of flat leaves of A. thaliana (Perez-Perez et al. 
2010). Nevertheless, information on gene networks that might 
control these regulators remains elusive. 

AS7 and AS2 are key regulators of the formation of flat 
symmetric leaves. AS1 and AS2 encode nuclear proteins and 
form a complex (referred to as AS1-AS2 in this report; Xu et al. 
2003, Yang et al. 2008, Luo et al. 2012). Mutations in these genes 
are associated with pleiotropic abnormalities in leaves along the 
three developmental axes (Redei and Hirono 1964, Tsukaya and 
Uchimiya 1997, Byrne et al. 2000, Ori et al. 2000, Semiarti et al. 
2001, Iwakawa et al. 2002, Iwakawa et al. 2007), suggesting that 
AS1-AS2 regulates multiple genes (Takahashi et al. 2008) that 
might be involved in leaf formation along these axes. The AS1- 
AS2 complex directly represses the transcription of BP and 
KNAT2 (Guo et al. 2008). Some of the pleiotropic abnormalities, 
including short leaves, of as2 and asl plants have been attrib- 
uted to ectopic expression of BP, KNAT2 and KNAT6 (Ikezaki 
et al. 2010), suggesting a role for these KNOX genes in the 
proximal-distal development of the leaf. Furthermore, tran- 
scripts levels of the ARF3/ETT, KAN2 and YAB5 genes are re- 
pressed downstream of AS7 and AS2 in shoot apices (Iwakawa 
et al. 2007, Takahashi et al. 2008). AS1-AS2 directly represses 
expression of ARF3/ETT, and the up-regulated ARF3/ETT in asl 
and as2 is responsible for defects in both development of the 
adaxial domain and expansion of the leaf lamina (Iwasaki et al. 
2013). These results suggest the involvement of AS1-AS2 in 
both adaxial development and the expansion of leaves through, 
at least in part, the functions of ARF3/ETT. 

Defects in polarity of asl and as2 leaves are enhanced under 
certain growth conditions as well as in conjunction with mu- 
tations in members of certain groups of genes (see the 
Introduction of Horiguchi et al. 2011b, Kojima et al. 2011, 
Ishibashi et al. 2012, Nakagawa et al. 2012, Xu et al. 2012), 
which are designated as modifiers of adaxial-abaxial patterning 
(Szakonyi et al. 2010, Iwasaki et al. 2013). These modifier genes 
include several that mediate the biogenesis of tasiR-ARF [a 
trans-acting short interfering RNA (ta-siRNA) see above]. 
Other relevant modifier genes belong to several different 
groups: those for ribosome biogenesis (Pinon et al. 2008, Yao 
et al. 2008, Horiguchi et al. 2011a); chromatin modification 
(Ueno et al. 2007, Kojima et al. 2011); genomic stability 
(Inagaki et al. 2009); and cell proliferation (Yuan et al. 2010, 
Horiguchi et al. 2011b, Ishibashi et al. 2012, Xu et al. 2012). We 
have previously shown that mutations of the gene encoding the 
histone acetyltransferase ELONGATA3 (EL03), a component of 
the Elongator complex, enhanced the polarity defects asso- 
ciated with as2 (Kojima et al. 2011). The as2 elo3 double 
mutant produces filamentous leaves with abaxialized epider- 
mis. We have further reported a mutation, designated enhancer 
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of asymmetric leaves! and asymmetric leaves! (eal), that is 
associated with efficient generation of abaxialized filamentous 
leaves on the as2 or asl background (Ishibashi et al. 2012). EAL 
encodes BOBBER1 (BOB1) (Jurkuta et al. 2009, Perez et al. 2009), 
an Arabidopsis ortholog of eukaryotic NudC domain proteins. 
Transcript levels of ARF3/ETT, ARF4, KAN2, YAB5 and all four 
class 1 KNOX genes are markedly elevated in shoot apices of asl 
eal and asl elo3 mutants. While these observations do suggest 
genetic interactions between AS1 and AS2 and each of these 
modifier genes, our understanding of the regulation mechan- 
ism for the expression of polarity-related effectors by AS1-AS2 
is still limited. Pathways of regulation by modifiers for the es- 
tablishment of leaf polarity and cell proliferation to produce flat 
and symmetric leaves are largely unknown. 

In the present study, we have carried out clustering analysis 
by KB-FuzzyART with a new gene list, Gene-list-2, by which we 
here report that AS1-AS2 plays roles in the repression of the 
Isopentenyltransferase (IPT) gene for cytokinin biosynthesis and 
that of the Kip-related protein (KRP) genes for inhibitors of 
cyclin-dependent kinases (CDKs) through the functions of 
ARF3/ETT and ARF4. In addition, our results suggest that 
KB-FuzzyART has high potential to extract candidate genes 
from among thousands of genes by using gene expression 
data of mutant strains. 



Results 




ARF3/ETT is a common target of modifier genes, 
BOB1/EAL and EL03/EAST1, in the asl 
background 

As mentioned in the Introduction, mutations of BOB1/EAL and 
EL03/EAST1 act as modifiers of asl and asl to enhance both 
the adaxial defect to generate filamentous leaves with an abax- 
ialized surface, and the expression levels of many genes includ- 
ing ARF3/ETT increased in asl bobl/eal and asl e\o3 double 
mutants (Kojima et al. 2011, Ishibashi et al. 2012). In addition, 
introduction of the ett mutation into the asl-1 eal~l double 
mutant partially suppressed its phenotype of filamentous 
leaves (Ishibashi et al. 2012), indicating that the increased 
level of ARF3/ETT transcripts is responsible for the adaxial 
defect. To confirm this indication further, we examined the 
effects of the ARF4 gene mutation, which is functionally redun- 
dant with ARF3/ETT (Pekker et al. 2005, Iwasaki et al. 2013), in 
the leaf phenotype of asl-1 eal-1 ett-13. As shown in Fig. 1Hand 
Table 1, the asl-1 eal-1 ett~13 arf4~l quadruple mutant plants 
generated symmetrically expanded leaves and no filamentous 
or trumpet-shaped leaves: phenotypes were more clearly res- 
cued in the quadruple mutant than in the asl-1 eal-1 ett~13 
triple mutant. 

We subsequently ascertained the effect of the ett-13 muta- 
tion on filamentous leaf formation in the asl-1 elo3-!7/east1-1 
double mutant, because the transcript level of ARF3/ETT was 
also increased in the double mutant (Kojima et al. 2011). As 
shown in Fig. 1K and Table 1, 70% of the asl-1 elo3-17 ett-13 



triple mutants exhibited expanded leaves. As shown in Fig. 1L 
and Table 1, the asl-1 elo3-U ett-13 arf4-1 quadruple mutant 
plants generated expanded leaves and no filamentous or 
trumpet-shaped leaves: the phenotypes were more clearly res- 
cued in the quadruple mutant than in the asl-1 elo3-17 ett-13 
triple mutant. Thus, the mutations of ARF3/ETT and ARF4 
clearly rescued filamentous leaves generated by asl-1 elo3-17, 
although the partner in double mutations was different 
from eal-1. The results of experiments with eal-1 and elo3-17 
suggest that the repression of ARE 3 /ETT and ARF4 is generally 
critical for establishment of the adaxial domain during leaf 
development. 

Meta-analysis to extract candidate downstream 
genes of AST, AS2 and ARF3/4 

We hypothesized that the population of up-regulated genes in 
asl-1 eal-1 and asl-1 elo3-17 might include candidate effector 
genes more closely related to leaf abnormality, because ARF3I 
ETT and ARF4 appear to encode transcription factors. Since 
levels of transcripts of abaxial determinant genes, such as 
ARF3/ETT and ARF4, are increased in asl-1 eal-1 and asl-1 
elo3-17 double mutants (Kojima et al. 2011, Ishibashi et al. 
2012), we can expect to extract efficiently candidate effector 
genes downstream of AS1, AS2, ARF3/ETT and ARF4 by clus- 
tering analysis of microarray data using these double mutants. 

We have already examined the gene expression profiles of 
Col-0, as1-1, asl-1 and ASl-eoe/ asl-1, in which AS2 was ecto- 
pically expressed (see the Materials and Methods) (Takahashi 
et al. 2008), by microarray analysis, and clustered genes accord- 
ing to their expression patterns using KB-FuzzyART, a clustering 
algorithm for the analysis of gene expression profiles (Takahashi 
et al. 2008). We re-calculated data sets of these plants by im- 
proved KB-FuzzyART described bellow. We named these data 
Dataset-A (Fig. 2A). In the present study, we examined the 
gene expression profiles of both asl-1 eal-1 and asl-1 elo3-17 
double mutants by microarray analysis, and clustered identified 
genes according to their expression patterns by KB-FuzzyART. 
We named these data Dataset-B (Fig. 2A). We then searched 
for candidate effector genes downstream of AS1, AS2, ARF3I 
ETT and ARF4 by meta-analysis using Dataset-A and Dataset-B 
(Fig. 2A). 

Clustering analysis by improved KB-FuzzyART 
using Gene~list~2 

In the previous report, we constructed and used Gene-list-1 as 
knowledge required for KB-FuzzyART. Since Gene-list-1 did not 
include some genes related to development and cell prolifer- 
ation, we added 103 genes to update Gene-list-1 and renamed 
it Gene-list-2 (Table 2, Supplementary Table S1) (see the 
Materials and Methods). We examined the gene expression 
profiles of Col-0, asl-1 eal-1 and asl-1 elo3-27 by microarray 
analysis, and clustered genes according to their expression pat- 
terns using KB-FuzzyART. As shown in Fig. 2B, we used data 
sets of Col-0, asl-1 eal-1 and asl-1 elo3-17 (Dataset-B). The 425 
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Col-0 (WT) 



as2-1 



ett-13 



arf4-1 




eal-1 



as2-1 
eal-1 



as2-1 
eal-1 
ett-13 



as2-1 
eal-1 
ett-13 arf4- 



* mi m * 



elo3-27 



as2-1 
elo3-27 



as2-1 
elo3-27 
ett-13 



as2-1 
elo3-27 
ett-13 arf4-1 



Fig. 1 The ett and ar{4 mutations efficiently suppressed the abaxial-adaxial defect phenotype of as2~l eal-1 and as2~1 elo3~27 leaves. Gross 
morphology of (A) Col-0, (B) aslA, (C) ett-13, (D) arf4-1, (E) eal-1, (F) as2-1 eal-1, (G) as2-1 eal-1 ett-13, (H) as2-1 eal-1 ett-13 arf4-1, (I) elo3-21, (J) 
as2-1 elo3-21, (K) as2-1 elo3-21 ett-13 and (L) as2-1 elo3-21 ett-13 arf4-1 plants is shown. Scale bars: 5 mm. Plants shown in A-L were photo- 
graphed 28 d (A-H) and 21 d (l-L) after sowing. White arrowheads and arrows indicate filamentous leaves and trumpet-shaped leaves, 
respectively. Higher magnification views of filamentous leaves and trumpet-shaped leaves are shown in F, J and K. Scale bars: 1 mm in higher 
magnification views. 



genes in Gene-list-2 were filtered and excluded 299 genes that 
have undetectable calls (Abs.), no-change calls (N.C.) or no 
probes on the ATH1 chip against the wild type in these 
plants. The expression data sets of the 126 filtered genes were 
processed by KB-FuzzyART and classified into 14 clusters, desig- 
nated as clusters B1-B14 (Fig. 2B). Then 22,746 probes on the 
GeneChip were filtered and 17,123 genes were excluded. The 



remaining 5,623 genes were assigned either as members of each 
cluster or as outliers. The patterns of gene expression levels in 
the outliers did not match those in any clusters. The outliers 
comprised 199 of the 5,623 genes (Table 3). Profiles of the 
constructed clusters and assignment of genes for Dataset-B 
are shown in Table 3 and Fig. 3A, B. To search for common 
genes downstream of EAL and EL03 on the 0s2-7 eal-1 and as2-1 
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Table 1 Frequencies of plants with filamentous and trumpet- 
shaped leaves 



Genotype 


No. of 
plants 
examined 


Filamentous 
leaves 


Trumpet- 
shaped 
leaves 


Col-0 


64 


0 (0%) 


0 (0%) 


as2-1 


91 


0 (0%) 


0 (0%) 


eal-1 


132 


0 (0%) 


0 (0%) 


as2-1 eal-1 


111 


84 (76%) 


16 (14%) 


as2-1 eal-1 ett-13 


45 


15 (33%) 


8 (18%) 


as2-1 eal-1 ett~13 arf4~1 


44 


0 (0%) 


0 (0%) 


elo3-27 


49 


0 (0%) 


0 (0%) 


as2-1 elo3-27 


44 


44 (100%) 


35 (80%) 


as2-1 elo3-27 ett-13 


44 


13 (30%) 


11 (25%) 


as2-1 elo3-27 ett-13 arf4-1 


5 


0 (0%) 


0 (0%) 



Frequency is defined by the ratio of the number of plants with more than one 
filamentous or trumpet-shaped leaf to the total number of plants examined. 
The percentages of plants with filamentous or trumpet-shaped leaves are indi- 
cated in parentheses. Plants were grown at 22°C. 

elo3-27 background, we focused on genes in the clusters B1, B6 
and B13, in which transcript levels of the genes were increased 
in both as2-1 eal-1 and as2A elo3-27 double mutants (Table 3, 
Fig. 3A, B). The cluster B1 contained 140 genes including ARF4. 
The cluster B6 contained 93 genes including BP/KNAT1, YAB5 
and KAN2. The cluster B13 contained 63 genes including ARF3/ 
ETT (Fig. 3, Table 3, Supplementary Table S2). 

We re-calculated data sets of Col-0, as 1 -1, as2-1 and 
AS2-eoe/as2-1 (Dataset-A) by KB-FuzzyART using Gene-list-2. 
The 425 genes in Gene-list-2 were filtered and the remaining 
114 genes were processed by KB-FuzzyART and classified into 
18 clusters, designated as clusters A1-A18. Then 22,746 probes 
on the GeneChip were filtered and the remaining 4,759 genes 
were assigned as members of each cluster or as outliers. The 
outliers comprised 62 of the 4,759 genes (Supplementary 
Table S3). Profiles of the constructed clusters and assignment 
of genes for Dataset-A are shown in Supplementary Table S3 
and Fig. S1A, B. 

It was previously reported that AS1 and AS2 form a repres- 
sor complex that binds directly to the 5' upstream regions of 
the KNOX loci BP/KNAT1 (Guo et al. 2008). We have also found 
that the AS1-AS2 complex binds to the 5' upstream regions of 
the ARF3/ETT locus and represses its expression (Iwasaki et al. 
2013). Therefore, AS1-AS2 is assumed to function as a tran- 
scriptional repressor and we focused on cluster A3, which 
included ARF3/ETT and BP/KNAT1. The expression levels of 
genes in cluster A3 were up-regulated in the as1-1 and 0s2-7 
plants, and down-regulated in the AS2-eoe/as2-1 plants 
(Supplementary Table S3). Cluster A3 contains 38 genes 
that include the BP/KNAT1, KNAT6, ARF3/ETT and YAB5 
genes (Supplementary Table S4). Expression levels of the 
genes included in cluster A6 were decreased in as1-1 and 
as2-7, but not changed in AS2-eoe/as2-1 (Supplementary 
Table S3), and contained 319 genes including KAN2 and 
ARF4 (Supplementary Table S4, Fig. S1A, B). Therefore, we 



focused on genes classified into cluster A6 in addition to cluster 
A3 in Dataset-A (Fig. 2 A). 

Search for candidate genes downstream of 
ARF3/ETT and ARF4 by using KB-FuzzyART 
and meta-analysis 

We combined the 357 genes in clusters A3 and A6 of Dataset-A 
and the 296 genes in clusters B1, B6 and B13 of Dataset-B, and 
extracted 57 genes (Fig. 2A, Table 4). These 57 genes are 
common downstream genes among AS7, AS2, BOB/EAL and 
EL03. In addition, the 57 genes include candidates of genes 
regulated by ARF3/ETT and ARF4. We found that only KRP2 
and KRP5 among the 57 genes have been known as genes that 
were involved in cell division, the transcript levels of which were 
commonly increased in as1, as2 and these double mutants. 
Since we focused on the cell division-related genes that were 
downstream of AS1 and AS2, we further confirmed the tran- 
script levels of the KRP2 and KRP5 genes in the double mutants 
by quantitative real-time RT-PCR (qRT-PCR). As shown in 
Fig. 4A, transcript levels of these genes increased by 2- to 
4-fold in 0S2-7 eal-1. Fig. 4B shows that levels of the KRP2 
and KRP5 gene transcripts increased by 1.5-fold in the 0s2-7 
elo3-27 mutant. Our results suggested that these two genes are 
negatively controlled downstream of AS2, EAL7 and EL03. 

We further examined effects of the ett-73 and arf4-1 muta- 
tions on increased transcript levels of these genes in the as2-1 
eal-1 mutant. As shown in Fig. 4C, transcripts levels of KRP2 
and KRP5 in the 0S2-7 eal-1 att-13 arf4-1 quadruple mutant 
decreased to levels close to those of both the wild-type 
plant and the att-13 arf4-1 mutant. These results suggest that 
KRP2 and KRP5 were positively regulated downstream of ARF3I 
ETT and ARF4 in the 0S2-7 eal-1 mutant. 

Since we have shown that shoot regeneration from leaf sec- 
tions in asl and as2 and the frequencies of shoot regeneration 
from leaf sections of the 0s7-7 and 0S2-7 mutants were restored 
by introduction of an ett-73 arf4-1 double mutation to that of 
the wild-type leaf sections, we speculated that some genes 
related to cytokinin synthesis and(or) signaling might be down- 
stream of AS1-AS2-ARF3/4. We did not find genes related to 
cytokinin, however, in 57 common candidate downstream 
genes of both the clusters A3 and A6 in Dataset-A and the 
clusters B1, B6 and B13 in Dataset-B. Therefore, we searched 
for genes related to cytokinin synthesis and (or) signaling in 
clusters B1, B6 and B13 and found the IPT3 gene in cluster B13. 
The IPT genes encode adenylate isopentenyltransferase, a cyto- 
kinin biosynthesis enzyme, in A. thaliana (Kakimoto 2001, Takei 
et al. 2001, Sakakibara et al. 2005). We further confirmed the 
transcript levels of IPT3 in the double mutants by qRT-PCR. As 
shown in Fig. 4A and B, transcript levels of IPT3 increased by 
3-fold in as2-1 eal-1 and by 1.5-fold in as2-l elo3-27. We also 
examined the effects of the ett-73 and arf4-l mutations on 
increased transcript levels of IPT3 in the as2-l eal-1 mutant. 
Fig. 4C shows that transcript levels of IPT3 in the as2-1 eal-1 
att-13 ar}4-1 quadruple mutant decreased. These results 
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Fig. 2 Gene expression analysis was conducted by the pre-processing 
procedure and clustering analysis by KB-FuzzyART. (A) RNA 
purification was conducted from the shoot apices of 15-day-old 
plants for Dataset-A (e.g. Col-0, as 1-1, as2-l and AS2-eoe/as2'1 mu- 
tants) and of 14- or 17-day-old plants for Dataset-B (e.g. Col-0, c?s2-7 
eal~l and as2-7 elo3~27 mutants). Gene expression data were 
obtained by using Affymetrix GeneChip ATH1. Each array data set 
comprised 22,746 A. thaliana genes (probe sets). KB-FuzzyART ana- 
lysis was conducted for two data sets, and gene clusters were con- 
structed. Furthermore, meta-analysis (integration of several analyses) 
was applied to clustering results for extraction of candidate 
downstream genes of AS 7 and AS2. (B) Clustering analysis of 
gene expression data (Dataset-B) was conducted by KB-FuzzyART. 
Each array data set and the 425 genes listed in Gene-list-2 were 
filtered by a pre-processing procedure (see the Materials and 
Methods), (i) First step: clustering of listed genes (knowledge-based 
listed genes) (e.g. cell fate determination genes, cell division-related 
genes, plant hormone-related genes), (ii) Second step: clustering of 
all genes was conducted by assignment to the 14 clusters constructed 
for Dataset-B. On the basis of this procedure, clustering by 
KB-FuzzyART provides the advantages of fast calculation and easy 
interpretation. 



Table 2 Categories and numbers of genes in Gene-list-2 



Category 


Subcategory 


No. of 






genes 


Cell fate 


Meristem-related genes 


17 


determination- 


Adaxial determinants 


6 


related genes 


Abaxial determinants 


11 




Organ differentiation-related genes 


19 




Small RNA biogenesis-related genes 


24 




AS2/LOB family 


41 




TCP family 


17 




Chromatin-related genes 


8 


Cell division 


Cell division-related genes 


61 


related genes 






Plant hormone- 


TIR1 gene and its homologs 


A 

4 


relatedgenes 


AUX/IAA family 


29 




ARC famik/ 
/Alxr Idllllly 


21 




GH3 family 


20 




SAUR family 


69 




IPT genes 


9 




AHK genes 


4 




ARR family 


24 




PIN genes 


8 




ACS family 


12 




AHP family 


6 




YUCCA family 


5 


Unclassified genes 


LSH family 


10 


Total 




425 



suggest that IPT3 was a positively regulated gene downstream 
of ARF3/ETT and ARF4 in the gs2-7 eal-1 mutant. Although we 
further examined transcript levels of /PT5 and IPT7, we 
observed no changes of transcript levels in either the 0S2-7 
eal-1 or the 0s2-7 elo3-27 mutant (data not shown). 

AS1-AS2 acts as the repressor of multiple 
downstream genes in independent pathways 

The BP/KNAT1, ARF3/ETT and YAB5 genes encode transcrip- 
tion factors, the target genes of which have yet to be identified. 
We here examined whether AS7 and AS2 might repress the 
expression of these genes in independent pathways, or whether 
some of them might control others downstream of AS7 and 
AS2. To test these possibilities, we investigated the effects of 
mutations in BP/KNAT1, KNAT2, KNAT6, ARF3/ETT, ARF4 and 
YAB5 on the transcript levels of BP/KNAT1, YAB5 and ARF3I 
ETT on the #s7 or as2 background by qRT-PCR (Fig. 5). As 
shown in Fig. 5A, transcript levels of BP/KNAT1 and YAB5 were 
higher in the 0S2-7 ett-13 arf4-l mutant than in the ett-13 arf4-l 
mutant. These results indicate that AS7 and AS2 repress BP/ 
KNAT1 and YAB5 independently from the ARF3/ETT and ARF4 
functions. 

We next measured levels of BP/KNAT1 and ARF3/ETT tran- 
scripts in the yabS-1, asl-1 yab5-l and as2A yab5-l mutants. 
The levels of BP/KNAT1 and ARF3/ETT transcripts were higher 
in the gs7-7 yab5-l and 0s2-7 yc?b5-7 mutants than in theyab5-7 
mutant (Fig. 5B). These results indicate that repression of BP/ 
KNAT1 and ARF3/ETT by AS 7 and AS2 is independent of YAB5 
functions. 
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Table 3 The constructed clusters and assignment of genes for 



Dataset-B 


Cluster 
number 


Cluster 
member 


Pattern 




as2-1 eal-1 


as2-1 elo3-27 


B1 


140 


+ 


+ 


B2 


459 


+ 


0 


B3 


1,173 


- 


- 


B4 


55 


- 


+ 


B5 


603 


0 


+ 


B6 


93 


+ 


+ 


B7 


552 


- 


0 


B8 


1,195 


+ 


0 


B9 


695 


0 




B10 


172 






B11 


43 


+ 




B12 


129 






B13 


63 


+ 


+ 


B14 


52 


+ 




Outlier 


199 






Overall 


5,623 







Cluster B1 (140 genes) 
(ARF4, KRP2) 



+, - and 0 indicate up-regulation, down-regulation and no change, respectively. 

Finally, we measured levels of ARF3/ETT and YABS tran- 
scripts in the bp- 7 knat2-3 knat6-2 and as2-l bp-1 knat2-3 
knatG-2 mutants. Transcript levels of ARF3/ETT and YAB5 
were higher in as2A bpA knat2-3 knat6-2 than in bp-1 
knat2-3 knat6-2 (Fig. 5C). These results indicate that AS7 and 
AS2 repress ARF3/ETT and YAB5 in a manner independent from 
the functions of BP/KNAT1, KNAT2 and KNAT6. 

Taken together, the results suggest that AS1-AS2 represses 
transcript levels of BP/KNAT1, ARF3/ETT and YABS genes in 
pathways independently of each other (Fig. 6). 



Discussion 

The AS1-AS2-ETT pathway controls gene 
networks for stabilization of the adaxial-abaxial 
patterning of leaves 

In the present study, we have shown that the phenotypes of 
filamentous leaves with adaxial defects that were generated by 
as2 bob1-3/eal-1 and as2 elo3-27 were rescued by introduction 
of the ett ar}4 double mutation to generate flat symmetric 
leaves (Fig. 1). We have previously shown that the phenotype 
of filamentous leaves of asl and as2 plants treated with 
5-aza-2 / -deoxycytidine, an inhibitor of DNA methylation, is 
completely suppressed by the introduction of ett-73 arj4A 
double mutations (Iwasaki et al. 2013). These results suggest 
that elevated expression of ARF3/ETT was responsible for the 
defects in adaxial development and the lack of lateral expan- 
sion. Therefore, repressive control of ARF3/ETT by AS1-AS2 
should be critical for stabilization of the normal development 
of leaves with proper adaxial-abaxial polarity and left-right 
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Fig. 3 Gene members were assigned by KB-FuzzyART clustering for 
Dataset-B. (A) Bar graph of genes as distributed within the 14 clusters 
of 5,623 genes. ARF4 and KRP2 were assigned to cluster B1. BP/KNAT1, 
YABS and KAN2 were assigned to cluster B6. ARF3/ETT, KRPS and IPT3 
were assigned to cluster B13. The blue portion indicates the number of 
transcription factors. (B) Clustering analysis by KB-FuzzyART. Levels of 
gene expression in comparison with those of Col-0 (WT) are shown. 
Clusters B1, B6 and B13 contained genes having elevated transcription 
levels in both c?s2-7 eal-1 and as2A elo3-27. 



symmetry. It might be possible that filamentous leaves gener- 
ated in many other kinds of modifier mutants (see Introduction 
of Kojima et al. 2011) in the asl or as2 genetic background 
might have resulted from elevated expression of both the 
ARF3/ETT and ARF4 genes. In consideration of these possibili- 
ties, we treated 0s2-7 ett- 7 3 arj4A triple mutant plants with a 
specific inhibitor of histone deacetylase (HDAC), trichostatin A 
(TSA), since 0S2-7 treated with TSA generates abaxialized fila- 
mentous leaves (Ueno et al. 2007). The phenotype of filament- 
ous leaves in 0s2-7 treated with TSA was effectively decreased in 
the ett-73 arf4-l mutant background (data not shown). 
Expression of ARF3/ETT is also regulated by several distinct 
mechanisms at the levels of transcription and translation 
(Nishimura et al. 2005, Ueno et al. 2007, Inagaki et al. 2009, 
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Table 4 The genes extracted commonly by using Dataset-A and 
Dataset-B 

AGI code Cluster Gene symbol Subcategory Cluster 

numbers name in numbers 

assigned Gene-list-2 assigned 

by using by using 

Dataset-B Dataset-A 



ATino73so 

AA I I vJU/ DJU 


B1 






A6 


ATin6309o 

AA I IVJUjUiU 


B1 


NRPD1A 

l\r\l U i AA 




A6 


AT1G64710 


B1 






A6 


AT1CI76R90 


B1 


GT2 




A6 


ATiri7R670 

AA I Ivj/OO/U 


B1 


AA / VjUnj 




A6 


AT1CI7RQ70 

AA I I Vj / 07/ U 


B1 


LUP1 




A6 


aticirociaO 

AA I IVJOUJJU 


B1 






A6 


AT9ri963nn 

AA I ZvjZOjUU 


B1 


CP Al PHA 7 

vj/ AAL/ FJAA / 




A6 


AT?ri?Q^^f) 

/A I £VJ£/JJU 


B1 






A6 


/A I I JU 


B1 






A6 


AT?r;4fic;Qn 

/A I iVJtDjyU 


B1 


DAG2 




A6 


AA I JVJ JUffU 


B1 


A/1 ES 70 




A6 


Aa I JUjUO JU 


B1 




Cell divisiorvrclcitGcl genes 


A6 


AT £ iCI01R70 

/A 1 JVJ VJ IO/U 


B1 






A3 


/A 1 JVJUZ.VJ.JlJ 


B1 


BLR/RPL 


/MprKf~pm-rpl;3l~prl crpnpQ 

JVlClloLCIII 1 CldLCU gtl ICo 


A6 


AA I JvJH/ IOU 


B1 






A3 


AA I JVJ JJ'tUU 


B1 


GLT1 




A6 


AA I JvJ J/ / I U 


B1 






A6 


AA I JVJOU'TjU 


B1 


ARF4 


Abcixicil determ i n3.nts 


A6 


AT5G61420 


B1 


MYB28 




A6 


AA 1 1 UZOJOU 


B6 






A6 


ATin 39940 

AA I IvJjZZ'-rU 


B6 


KAN 2 


A Kc» y id 1 Hptprm inantc 

AAUd.XId.1 Ut.Lt: I 1 1 III Idl lLo 


A6 


aa i i vj jj / y\j 


B6 






A6 


AT1 CI 77760 

AA 1 1 VJ / / / OU 


B6 


NIA1 




A6 


AT? CI f) 1670 


B6 


ntnuHtl 7 




A6 


AT9n99660 

AA 1 ZvJZZDDU 


B6 






A6 


AT9CI96c;R0 

AA 1 ZVJZUJOU 


B6 


YA B5 


A Kd Yid 1 Hpfprm inanfc 

/AL/dAldl UtLtl 1 I III 1 dl ILo 


A3 


AT9CI9Q660 

AA 1 Z. Vj ZyUUU 


B6 






A6 


AT9n4ACi10 

AA I iVJ'iJJ IU 


B6 


ATT/ 7 




A6 


ATaCI 09390 

AA I JvJU7J7U 


B6 


MT2A 




A6 


AT^ri9^9Q0 

AA I JUZjZyU 


B6 


/_SH4 


LJM Idl 1 Illy 


A6 


ATaCIc;9R00 

AA I JVJ JiOUU 


B6 






A6 


AT4CI09 3R0 

AA I fVJU^JOU 


B6 


SAG27 




A6 


AT4G08150 


B6 


BP/KNAT1 


Meristem-related genes 


A3 


AT4G 16447 


B6 






A6 


AT4G34710 


B6 


ADC2 




A6 


AT4G38210 


B6 


ATEXPA20 




A6 


AT4G 38620 


B6 


MYB4 




A6 


AT5G 14090 


B6 






A6 


AT5G39850 


B6 






A6 


AT5G43580 


B6 






A6 


AT5G64110 


B6 






A6 


AT1G05710 


B13 






A6 


AT2G05380 


B13 


GRP3S 




A6 


AT2G33860 


B13 


ARF3/ETT 


Abaxial determinants 


A3 


AT2G 34650 


B13 


PID 




A6 


AT2G 38870 


B13 






A6 



(continued) 



Table 4 Continued 



AA VJ 1 LUUC 


numbers 
assigned 
by using 
Dataset-B 


{"■pi"ip c\#m Hnl 
vjciic syiiiuui 


juuLaicgui y 

name in 
Gene-list-2 


numbers 
assigned 
by using 
Dataset-A 


AT3G02500 


B13 




- 


A3 


AT3G 13790 


B13 


ATBFRUCT1 


- 


A6 


AT3G22550 


B13 




- 


A6 


AT3G23050 


B13 


AXR2/IAA7 


AUX/IAA family 


A6 


AT3G24810 


B13 


KRP5/ICK3 


Cell division-related genes 


A6 


AT4G00720 


B13 


ATSK32 




A6 


AT4G32150 


B13 


ATVAMP711 




A6 


AT4G 37410 


B13 


CYP81F4 




A6 


AT5G40540 


B13 






A6 


AT5G65310 


B13 


ATHB5 




A6 



-, the gene was not included as a member of Gene-list-2. 



Ng et al. 2009, Iwasaki et al. 2013). In addition, mRNAs of ARF3/ 
ETT and ARF4 are both targeted by a ta-siRNA called tasiR-ARF 
(Allen et al. 2005, Williams et al. 2005, Fahlgren et al. 2006, 
Hunter et al. 2006, Nogueira et al. 2006, Nogueira et al. 2007, 
Chitwood et al. 2009, Schwab et al. 2009). The presence of these 
multiple processes for ARF3/ETT repression implies that the 
negative control of ARF3/ETT (probably also ARF4) as well as 
class 1 KNOX genes is essential for the normal formation of 
vegetative and floral lateral organs in A. thaliana (Tabata 
et al. 2010). 

AS1-AS2 regulates KRP and IPT genes through 
ARF3/ETT and ARF4 for leaf development 

The present study has also shown that expression of KRP2, 
KRPS and IPT3 is positively controlled downstream of ARF3I 
ETT, at least on the as2 and modifier backgrounds, suggesting 
that ARF3/ETT has abilities to elevate the endogenous level of 
cytokinin and to repress cell cycle progression (Fig. 6). 

What is the developmental relevance of the repressive con- 
trols of these three genes by AS1-AS2? We postulate that re- 
pression of KRP2, KRPS and IPT3 by AS1-AS2 through ARF3I 
ETT functions might be related to the stabilization of leaf 
development (Fig. 6). 

KRP2 and KRPS of A. thaliana encode cyclin-dependent 
kinase inhibitors (CKIs) (De Veylder et al. 2001). Leaf formation 
might also require that cell proliferation that might be achieved 
by the properly repressive control of CKI genes by AS1-AS2 
(Fig. 6). Some CKI genes are predicted to be involved in endor- 
eduplication in A. thaliana, which eventually leads to cell cycle 
arrest (Ormenese et al. 2004). KRP2, for example, is expressed in 
tissues in which endoreduplication occurs, and its overexpres- 
sion (or ectopic expression) results in an increase in the plolidy 
level (Verkest et al. 2005). Our previous observation that the 
rate of endoreduplication is stimulated during leaf formation in 
as2A eal-1 double mutants (Ishibashi et al. 2012) is consistent 
with the results in the present study. There might be a correl- 
ation between the progression of endoreduplication and cell 
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Fig. 4 Levels of transcripts of KRP2, KRPS and /PT3 genes in shoot apices of various mutants measured by quantitative real-time RT-PCR (qRT- 
PCR). Transcript levels are shown for (A) Col-0, aslA, eal-1 and aslA eal-1 mutants; (B) Col-0, aslA, elo3-27 and as2A elo3-27 mutants; and (C) 
Col-0, as2-l eal-1, ett-13 arf4-l and as2-l eal-1 ett-13 arf4-l mutants. Total RNA of plants in (A) and (C) was extracted from short apices of 
seedlings 17 d after sowing. Total RNA of plants in (B) was extracted from short apices of seedlings 14 d after sowing. Each value was normalized 
by reference to the level of ACTIN2 transcripts. The values from Col-0 plants were set arbitrarily at 1.0. Bars indicate the SD among more than 
three biological replicates. 



differentiation in plants (Joubes and Chevalier 2000, De Veylder 
et al. 2002, Inze and De Veylder 2006). KRP2 is expressed in 
mitotically dividing cells (Ormenese et al. 2004). Therefore, 
AS1-AS2 might be involved in the modulation of a balance 
between cell cycle progression and cell differentiation by the 
proper repression of expression of these CKI genes. 

CKIs, which interact with CDKs to inhibit their kinase activ- 
ities, act as key regulators of cell cycle progression (Sherr and 
Roberts 1999). Typical examples are p21, p27 and p57 in animal 
cells (Toyoshima and Hunter 1994, Lee et al. 1995). A role for 
the CKIs in cell differentiation was reported during muscle de- 
velopment and maintenance of adult hematopoietic stem cells 
(Zhang et al. 1999, Matsumoto et al. 2011); however, little is 
known about how the cell cycle machinery communicates with 
developmental processes in both plants and animals. 
Understanding of molecular networks that involve CKIs must 
shed light on the coordination between cell proliferation and 
the developmental processes in multicellular organisms. 

The endogenous level of cytokinin might increase around 
the SAM in asl and as2 mutants because of the ectopic expres- 
sion of /PT3, which might affect developmental states of cells in 
the leaf primordia of asl and as2 mutants. In fact, the regen- 
eration potential of shoots in leaves of the asl and as2 mutants 



increases during in vitro culturing (Semiarti et al. 2001). This 
phenotype is consistent with the above prediction of an in- 
crease in the cytokinin level in these mutant leaves and has 
recently been shown to contribute to the elevated expression of 
ARF3/ETT in asl and as2 (Iwasaki et al. 2013). Thus, an unneces- 
sary increase of cytokinin should disturb the stabilization of the 
developmental status of leaves. For further understanding of 
developmental roles of the AS1-AS2-ETT-IPT3 pathway in 
cytokinin synthesis, endogenous levels of cytokinin in asl and 
as2 must be examined by further experiments. 

Materials and Methods 
Plant materials and growth conditions 

Arabidopsis thaliana ecotypes Col-0 (CS1092), ash! (CS3374) 
and as2-7 (CS3117) were obtained from the Arabidopsis 
Biological Resource Center (ABRC). We outcrossed 0s2-7 with 
Col-0 three times and asl-1 with Col-0 once, and used the 
progeny for our experiments (Semiarti et al. 2001). For analysis 
of phenotypes, seeds were sown in soil or on plates of MS 
(Murashige and Skoog) medium. After 2 or 3 d at 4°C in dark- 
ness, plants were transferred to a daily regime of 8 h of darkness 
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Fig. 5 Transcript levels of BP/KNAT1, YAB5 and ARF3/ETT genes in the shoot apices of various 15-day-old mutants as measured by quantitative 
real-time RT-PCR (qRT-PCR). Shown are the transcript levels of (A) BP/KNAT1 and YAB5 in ett-13 arf4~l, as1-1 ett-13 arf4~l and as2-l ett-13 
arf4-1 mutants; (B) BP/KNAT1 and ARF3/ETT in the yab5-1 er-105, as1~l yab5-1 er-105 and as2-l yab5-1 er-705 mutants; and (C) ARF3/ETT and 
YAB5 in the bp~l knat2~3 knat6~2, asl~l bp~1 knat2~3 knat6~2 and as2~l bp~1 knat2~3 knat6~2 mutants. Each value was normalized by reference to 
the level of ACTIN2 transcripts. The values from Col-0 plants were set arbitrarily at 1.0. Bars indicate the SD among more than three biological 
replicates. 
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Fig. 6 Model for regulation of ARF3/ETT and ARF4 by the AS1-AS2 complex in the early stage of leaf primordia in A. thaliana. Temporal 
regulation of the gene expression for establishment of leaf adaxial-abaxial polarity is shown in A. thaliana. The leaf primordium in the early stage 
has begun to grow outward and then acquired polarity in the radial dimension (Lynn et al. 1999). In the leaf primordium AS1-AS2 represses class 
1 KNOX and abaxial determinant genes, together with EAL and EL03. AS1-AS2 and modifiers regulate KRP and IPT genes through ARF3/ETT and 
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and 16 h of white light at 50 umol m _2 s _1 at 22°C, as described 
previously (Semiarti et al. 2001). Ages of plants are given in 
terms of numbers of days after sowing. The eastl-1 /elo3-27 
and bobl-3/eal-1 mutants were previously described (Kojima 
et al. 2011, Ishibashi et al. 2012). The seeds of ett- 13 and arf4-l 
were kindly provided by Dr. Yuval Eshed (Weizmann Institute 
of Science, Israel). The ett-13 allele contains a loss-of-function 
mutation with a T-DNA insertion at 18 bp downstream of the 
translation start site of ETT (Pekker et al. 2005). The arf-4-1 
allele is also a loss-of-function mutation with a Ds insertion 
(Parinov et al. 1999). Both ett-13 and arf4-l are presumptive 
null, or close to null, alleles. The yab5-l mutant, which has a 
stop codon mutation in the YABBY domain, was obtained from 
the Arabidopsis TILLING Project (Till et al. 2003). This plant has 
not only yab5-l but also er-705 mutations. These two muta- 
tions could not be segregated by backcrossing three times be- 
cause of the short distance between YAB5 and ER (<100kbp). 
The bp-1, knatl-3 and knat6-2 mutants were previously 
described (Ikezaki et al. 2010). The transgenic plant AS2-eoe/ 
as2-l, in which AS2 was ectopically expressed under the control 
of the AS 7 promoter (pAS1:AS2; designated AS2-eoe), was also 
previously described (Iwakawa et al. 2002, Iwakawa et al. 2007). 

Microarray analysis 

Arrays of Dataset-A were described previously (Takahashi et al. 
2008). For Dataset-B, shoot apices of Arabidopsis plants were 
harvested 14 or 17 d after sowing, and immediately frozen in 
liquid nitrogen and stored at -80°C. Total RNA extracted from 
each sample of shoot apices containing leaf primordia and 
young leaves was reverse transcribed, yielding double-stranded 
cDNA, which was transcribed in vitro in the presence of 
biotin-labeled nucleotides using a GeneChip 3'IVT Express Kit 
(Affymetrix Inc.), and purified. Labeled amplified antisense RNA 
was fragmented and hybridized to Affymetrix ATH1 GeneChip 
arrays for 16 h at 45°C according to Affymetrix protocols. 
Arrays were then washed on an Affymetrix Fluidics Station 
450 and measured for fluorescence intensity with an 
Affymetrix GeneChip Scanner 3000 7G. The microarray data 
of Dataset-A presented in this paper are available from the 
Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/ 
query/acc.cgi?acc=GSE44028) under the accession number 
GSE44028. 

Data processing 

Raw data processing was performed by using Affymetrix Gene 
Chip Operating Software (GCOS) (Version 1.4.0.036). We used 
two sets of gene expression data, Dataset-A obtained from 
strains Col-0, asl-1, as2A and /KS2-eoelas2-1, and Dataset-B 
obtained from strains Col-0, as2A eal-1 and 0s2-7 elo3-27. 
Each set of array data comprised 22,746 plant genes (probe 
sets). Initially, we calculated the expression signals for all strains 
and the log2(ratio) for three strains (gs7-7, c?s2-7 and AS2-eoe/ 
as2~l) and two strains (as2-l eal-1 and c?s2-7 elo3-21) against 
Col-0 by GCOS. In this experiment, we excluded 64 controls and 



2,177 genes subject to cross-hybridization, according to NetAffx 
Annotation (www.affymetrix.com). Furthermore, for each 
group set, we excluded those genes for which all strain data 
sets showed an absent call (i.e. the detection call determined by 
GCOS based on the P-value of the one-sided Wilcoxon 
signed-rank test; an absent call means P > 0.065, which is the 
default threshold in GCOS), because it indicates that the ex- 
pression signal was undetectable. We also excluded those genes 
for which all strain data sets showed a no change call [i.e. no 
change call means (1-0.006) > P > 0.006 for one-sided 
Wilcoxon signed-rank test], because no change indicates that 
the expression signal is almost equal to that of Col-0. Thus, 
4,759 and 5,623 genes were selected for Dataset-A and 
Dataset-B, respectively. Among these data, the log2(ratio) 
values that were >2-fold or <0.5-fold were rounded to 2-fold 
and 0.5-fold, respectively; the log2(ratio) values with no change 
were rounded to 0, to avoid category proliferation in clustering. 

Gene-list-2 for Arabidopsis thaliana 

In the present study, we used KB-FuzzyART. This is a powerful 
algorithm for clustering of biological data (Takahashi et al. 
2008); however, KB-FuzzyART is a knowledge-based system 
that requires some knowledge data. Therefore, we constructed 
Gene-list-1 for A. thaliana genes in a previous study (Takahashi 
et al. 2008). We added 103 genes (e.g. chromatin-related genes, 
ACS family, LSH family and PIN genes) to update Gene-list-1, 
and renamed it Gene-list-2, which comprises 425 genes 
(Table 2). Further detailed information of Gene-list-2 is 
shown in Supplementary Table S1. 

The clustering by KB-FuzzyART 

We applied KB-FuzzyART for two data sets: Dataset-A (Col-0, 
0S7-7, as2A and AS2-eoe/as2-1) and Dataset-B (Col-0, 0S2-7 
eal-1 and as2-l elo3-27). First, for KB-FuzzyART, FuzzyART 
was applied to any knowledge, such as the 114 and 126 genes 
obtained from Gene-list-2 for each data set, respectively. The 
KB-FuzzyART parameters (i.e. vigilance parameter and learning 
rate) were optimized by using an optimal clustering index (OCI) 
(Takahashi et al. 2003), and then 18 and 14 clusters were con- 
structed, respectively, for each data set in the present study. 
Weight vectors (representative patterns for each cluster) for 
the 18 and 14 clusters were constructed by FuzzyART, and 
the learning rate was set to 0, to fix the constructed weight 
vectors. FuzzyART with 18 or 14 fixed weight vectors was then 
applied to the 4,759 or 5,623 filtered genes for each data set. All 
of the 4,759 or 5,623 genes were assigned to each cluster or 
outliers (62 or 199 genes) for each data set. The outliers repre- 
sent genes with patterns that did not match patterns from 
Gene-list-2. 

Real-time RT-PCR 

Leaves and shoot apices of wild-type and mutant plants were 
harvested at 14 or 17 d after sowing, and immediately frozen in 
liquid nitrogen and stored at — 80°C Total RNA was isolated by 
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using the RNeasy Plant Mini Kit (Qiagen). For the analysis of 
RNA levels by real-time PCR, we prepared 5 ug of total RNA 
to convert to cDNA as described by Iwakawa et al. (2007). PCR 
was performed in the presence of the double-stranded 
DNA-specific dye Power SYBR Green (Applied Biosystems) 
with primers as described (Iwakawa et al. 2007, Takahashi 
et al. 2010). Amplification was monitored in real time by 
using the Applied Biosystems StepOnePlus Real-Time PCR 
system, according to the supplier's recommendations. The 
mean value of three technical replicates was normalized by 
the value of ACTIN2 transcripts. Primers used are listed in 
Supplementary Table SS. 



Supplementary data 

Supplementary data are available at PCP online. 
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